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Abstract 



Using data from the HEGRA air shower array, taken in the period from April 1998 
to March 2000, upper limits on the ratio I^/Icr of the diffuse photon flux J 7 to 
the hadronic cosmic ray flux Iqr are determined for the energy region 20 TeV 
to 100 TeV. The analysis uses a gamma-hadron discrimination which is based on 
differences in the development of photon- and hadron-induced air showers after the 
shower maximum. A method which is sensitive only to the non-isotropic component 
of the diffuse photon flux yields an upper limit of 1^/Iqr (at 54 TeV) < 2.0 x 10~ 3 
(at the 90% confidence level) for a sky region near the inner galaxy (20° < galactic 
longitude < 60° and | galactic latitude | < 5°). A method which is sensitive to 
both the isotropic and the non-isotropic component yields global upper limits of 
I-j/Icr (at 31 TeV) < 1.2 x 10~ 2 and Ij/I C r (at 53 TeV) < 1.4 x 10~ 2 (at the 
90% confidence level). 

Key words: Cosmic rays, Diffuse gamma-ray Galactic emission; Gamma rays: 
observations. 



1 Introduction 



The flux of the diffuse photon background radiation, its spectrum and its evo- 
lution with time are very interesting topics of astrophysics. In nearly all energy 
regions both an extragalactic and a galactic component of the diffuse radia- 
tion can be identified. The former component, which is essentially isotropic, 
gives information about processes and developments at large distances from 
our galaxy, ranging partly back to the early universe. Details about our own 
galaxy like the interstellar matter, radiation and magnetic fields as well as 
about the origin and the propagation of galactic cosmic rays can be deduced 
from investigations of the galactic diffuse emission. 

Direct measurements of the diffuse emission in the 30 MeV to 50 GeV energy 
range are provided by the EGRET experiment [1,2]. At high galactic lati- 
tudes an extragalactic component has been observed, a significant fraction of 
which is attributed to the direct emission from AGNs [3]. The main part of 
the diffuse emission is concentrated at the galactic disc and the inner galaxy, 
which therefore points to a galactic origin. The EGRET data below 1 GeV 
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can be well explained by interactions of cosmic-ray electrons and protons with 
the interstellar radiation field and with the interstellar matter [4]. The domi- 
nant interaction processes are n production by nucleon-nucleon interactions, 
inverse Compton (IC) scattering with low-energy photons and high-energy 
electron bremsstrahlung. Above 1 GeV the observed photon flux exceeds the 
model predictions [1]. A possible explanation for this discrepancy is that the 
average interstellar proton spectrum is harder than the spectrum observed in 
the local neighbourhood, leading to a harder 7r° and thus also to a harder 
diffuse gamma-ray spectrum [5-8]. The discrepancy may also be due to an 
underestimation of the inverse-Compton contribution, which starts to become 
relevant above 1 GeV [9-15]. The favoured hypothesis for an enhanced IC con- 
tribution is an average electron spectrum in the galaxy which is substantially 
harder than that measured locally [10,14,15]. An alternative or additional 
contribution could be gamma-radiation from supernova remnants and pulsars 
[16]. 

It is interesting to note that in models with increased inverse-Compton con- 
tribution [13,14] not only the excess of photons above 1 GeV can be explained 
but also the longitude and latitude profiles up to the galactic poles can be 
reproduced. This could be relevant for the experimental determinations of the 
isotropic extragalactic emission. 

In these models the IC contribution to the diffuse photon background domi- 
nates above 1 GeV. Depending on the cutoff energy for the electron injection 
spectrum a maximum of the differential photon flux (multiplied by E 3 ) is pre- 
dicted between 30 TeV and 200 TeV [1 1] . For the inner galaxy the flux at the 
maximum corresponds to a fraction of ~ 5 x 10~ 4 of the cosmic ray flux (see 
Fig. 11). The IC contribution is rather broad in galactic latitude and extends 
up to the galactic poles, and even there it is comparable to or greater than 
the expected isotropic extragalactic contribution [13,14]. 

The broad distribution in galactic latitude is due to the broad distribution 
of the interstellar radiation field, which is the target field for the IC process. 
In contrast, a model with enhanced 7r° production would predict a gamma 
radiation which is more confined to the galactic disc, because of the narrow 
distribution of the interstellar matter, the target for the 7r° production process. 
This difference may allow one to distinguish between the two hypotheses. 

Contributions to the extragalactic diffuse emission below 100 TeV are expected 
from the cascading of ultra-high-energy photons and electrons on the Cosmic 
Microwave Background (CMB), creating eventually a pile-up of gamma-rays in 
the 10 - 100 TeV region [17,18]. The ultra-high-energy photons and electrons 
may be decay products of supermassive particles [19], themselves radiated 
during collapse or annihilation of topological defects [20]. Alternatively they 
may result from the interaction of ultra-high-energy hadronic cosmic rays with 
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the CMB [21]. The expected level of gamma flux contributed by the cascading 
process is estimated to be 1CT 5 of the cosmic ray flux [21-23]. 

In the 1 TeV to 1 PeV energy range there are two experiments which claim 
the observation of a diffuse photon signal: [24] and [25]. The result reported 
by [24] is only a 2.8 sigma effect and may therefore well be considered as an 
upper limit. Moreover, there may be an additional uncertainty of the result 
due to uncertainties on the muon content of photon showers [26]. The gamma- 
ray excess reported by [25] is a 3.8 sigma effect, which may be correlated with 
the Loop 1 region. In a later publication by the same collaboration [27] only 
an upper limit of the diffuse gamma-ray flux is given. All other experiments 
in the 1 TeV to 1 PeV energy region give upper limits on the diffuse photon 
flux (see Tables 2 and 3 and Figs. 11 and 12). 

In the present analysis [28] an independent new determination of an upper 
limit on the diffuse photon background is performed in the 20 TeV - 100 TeV 
region, where only two other experiments have reported results [29,30], in 
addition to previous measurements with the HEGRA array [31,32]. 

The energy region below and around 100 TeV is of particular interest because 
on the one hand the extragalactic component of the diffuse photon background 
is expected to contribute predominantly in this energy region. On the other 
hand, according to some models [11], the ratio L,/Icr of the differential diffuse 
photon flux 7 7 of galactic origin to the cosmic ray flux Iqr has a maximum in 
the 100 TeV region (see Figs. 11 and 12). 

The paper is structured as follows: Section 2 describes the experimental details, 
Section 3 the Monte Carlo simulation. In Section 4 the definition of the data 
sample to be used in the analysis is given. The measurement of the upper 
limit on the diffuse photon flux is presented in Section 5 and the results are 
discussed in Section 6. 



2 Experimental Details 

2.1 The Detector 

The data were taken with the HEGRA air shower array located at 2200 m 
a.s.l. on the Canary island La Palma. Since April 1998 the array consisted 
of 97 wide-angle Cherenkov detectors (AIROBICC) and of 182 scintillation 
detectors. A fire in October 1997 had destroyed 39 AIROBICC detectors, all 
of which were rebuilt, and 65 scintillation detectors, of which only 4 were 
rebuilt. The two types of detectors measure the Cherenkov radiation and the 
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secondary particles (mainly photons, converted in a 5 mm thick lead plate 
above the 4 cm thick plastic scintillator sheet, and electrons) respectively, 
which are produced in air showers induced by primary cosmic rays (including 
gammas) in the atmosphere. The AIROBICC detectors were distributed over 
an area of 200 m x 200 m. This area contained completely the 150 m x 200 m 
area covered by the scintillator array. The detectors are described in detail in 
[33] and [34]. 

The experimental trigger required a time coincidence within 150 ns of > 6 
Cherenkov detectors (AIROBICC trigger) or of > 14 scintillation detectors 
(scintillator trigger). The trigger rate was 28 Hz on the average. 

The data used in this analysis were taken in the period from April 1998 to 
March 2000. Only data with an AIROBICC trigger have been considered. The 
total effective observation time amounts to ~ 1464 hours. 

2.2 Detector Calibration and Shower Reconstruction 

The detectors register the arrival time {time data) and the light flux {amplitude 
data) hitting the photocathode of the photomultipliers. For each detector the 
amplitude is given in terms of the number of ADC channels. In the case of the 
scintillator array the relative calibration of amplitudes between the individual 
detectors is done on the basis of the so-called '1-MIP peak'. The 1-MIP peak 
is the position of the maximum in the ADC spectrum corresponding to the 
detector response for 1 minimum ionizing particle [35]. For the Cherenkov 
array a relative calibration between the individual detectors is achieved by 
normalizing the high-amplitude tails of the individual ADC spectra to each 
other [35]. The absolute calibrations of the amplitudes are described in Section 
3.2. 

By the absolute calibrations the amplitudes are transformed into a density of 
Cherenkov photons impinging on the surface of the Cherenkov detectors or a 
density of secondary electrons and photons impinging on the surface of the 
scintillation detectors respectively. 

For each shower the shower direction, the position of the shower core and the 
Cherenkov light density p{r) as a function of the distance r from the shower 
core were reconstructed by a simultaneous fit {global AIROBICC fit) to 
the time and amplitude data from the Cherenkov detectors. From p(r) the 
Cherenkov light density L 9Q at r = 90 m and the light radius Rl = 
— 1 / {dlnp{r)/dr) in the region 50 m < r < 120 m were derived. The total 
number N s of photons and electrons at the detector level was determined 
by a maximum likelihood fit to the particle densities from the scintillation 
detectors. 
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L90 is known to be a good measure of the shower energy. In the following 
L 90 and Rl are assumed to be given in units of [no. of photons / m 2 ] and [m] 
respectively. 

The quality of the various fits to the time and amplitude data, expressed in 
terms of x 2 — probabilities and in terms of errors of fitted parameters, is used 
below in the definition of the selection criteria for the showers. 



3 Monte Carlo Simulation and Data-Monte Carlo Comparison 

3.1 Monte Carlo Simulation 



An essential part of the present analysis is the Monte Carlo (MC) simulation 
of the experiment. The simulation comprises : 

- The simulation of the physical processes in the shower development, includ- 
ing the production of Cherenkov photons. The CORSIKA program, version 
4.068, was used for this step [37]. 

- The simulation of the detector. A comprehensive discription is given in [38]. 

- The simulation of the chemical composition and of the energy spectrum of 
the charged cosmic rays. The simulation was based on the compilation [39] 
(see Table 1). 

MC data were generated for primary photons, protons, He-, O- and Fe nuclei, 
for the zenith angles = 10°, 20° and 30° and for the energy region 5 TeV < 
E < 1000 TeV. For each particle type and each zenith angle approximately the 
same number of showers was generated, giving a total of 133 000 showers. By 
assigning to each of these showers 20 different core positions, chosen randomly 
in a 300 m x 300 m area centered at the actual detector array, the statistics 
of quasi-independent showers was increased by a factor of 20. 

The generated showers were given appropriate weights to simulate the ex- 
pected distributions in E, 9 and the cosmic ray species. For the comparisons 
with the experimental data an admixture of photons with O = 0.0002582 (m 2 - 
s-sr-TeV)" 1 has been assumed in the MC sample (see Table 1), corresponding 
to 1/1000 of the hadronic cosmic ray flux at 1 TeV. 

In the following 'MC events' always means 'weighted MC events'. 
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3.2 Absolute Calibration of N s and L 90 

The absolute calibration of N s and L 90 is done separately for each subrun 
(corresponding to 30 minutes of observation time). The calibration factor for 
N s is determined such that (log 10 (iV s )) in the experimental data agrees with 
(log 10 (7V s )) in the MC data. This is done using samples defined by the standard 
selection criteria (see Section 4). 

The absolute N s calibration defines - through the MC simulation - the absolute 
energy scale in the experimental data. By comparing with other methods of 
calibrating N s , for example using the 1-MIP peak not only for the relative but 
also for the absolute calibration of N s [35], the error in the absolute energy 
scale is estimated to be in the order 20%. 

After calibration, the log 10 (iV s ) distributions of the experimental data and of 
the MC data agree fairly well (see Fig. 1). There is also reasonable agreement 
in the total number of showers: for a small subset of runs the total number of 
showers in the experimental data was compared with the corresponding MC 
number expected on the basis of the measured on-time and the all-particle 
cosmic ray flux from [39]. There is agreement within 7%. This can be regarded 
as a good consistency between experimental data and MC data, given the fact 
that already a miscalibration of N s by 10% would yield a discrepancy in the 
number of reconstructed showers of (1.10 17 — 1.0) = 17.6%. 

The absolute calibration of Lgo is done using the calibrated N s values, namely 
by determining a calibration factor for L go such that (log 10 (A r s /L 90 )) in the 
experimental data coincides with (log 10 (A^ s /L 90 )) in the MC data. The cal- 
ibration factor is determined separately for each of five bins in the variable 

3. 3 Adjustment of Errors in the Monte Carlo Simulation 

A comparison of the fluctuations of arrival times and amplitudes between 
experimental data and MC data gave good agreement for the data from the 
scintillation detectors. For the Cherenkov detectors the fluctuations were found 
to be significantly lower in the MC data than in the experimental data. The 
difference is explained on the one hand by an underestimation of the average 
light of the night sky background in the MC simulation. On the other hand, 
the difference may be due to subtle effects in the conversion chain from the 
photons to the final electronic signal, such as variations in the Winston cone 
reflectivity, non-uniformity in the quantum efficiency and the photoelectron 
collection efficiency of the photomultiplyers, etc. . Instead of repeating the 
MC simulation with a more realistic night sky background and with a refined 



7 



simulation of the detector non-uniformities, the fluctuations of arrival times 
and amplitudes from the Cherenkov detectors were randomly increased in 
the MC-generated data to agree with those in the experimental data. On the 
average the RMS of the Cherenkov photon densities and that of the Cherenkov 
photon arrival times was increased by a factor of 1.8. 

3.4 Data-MC Comparison 

Some comparison of experimental data and MC data are shown in Figs. 1 to 
4. The data samples are defined by the standard selection criteria listed in 
Section 4. Because these selection criteria are tighter than the trigger condi- 
tion no simulation of the trigger was necessary. In each of the Figs. 1 to 3 
the two distributions were normalized to the same total number of events. For 
log 10 (7Vg), Rl and also for log 10 (iV s /L 9 o) versus (— 1/Rl) the agreement be- 
tween data and MC is quite good. A large discrepancy is found for log 10 (L 90 ) 
at log 10 (L 90 ) < 4.0. This may be due to the time dependence of the light of 
the night sky in the experimental data which was not taken into account in 
the MC simulation. 

3.5 Comparison of photon-induced to hadron-induced showers in the MC data 

An experimental upper limit of the flux of diffuse photons will obviously be the 
tighter the better photons can be discriminated from hadrons. It is therefore 
worthwhile searching for measurable quantities which allow a photon-hadron 
discrimination at least to a certain degree. In the present analysis this is done 
exclusively with MC data: As no prominent photon signal has been found so 
far with HEGRA-array data there is no experimental data sample, sufficiently 
enriched with photons, that could be used for studying the photon-hadron 
separation. 

As shown in [41] useful variables in this context are 1/Rl and the ratio 
log 10 (A" s /L 90 ): because hadron-induced showers in general develop more slowly 
after the shower maximum than photon-induced showers, the ratio log 10 (A" s /L 90 ) 
at fixed position of the shower maximum (estimated by 1/Rl [40]) is expected 
to be on the average larger for hadron than for photon showers. This be- 
haviour is illustrated in Fig. 5. The photon-hadron difference is quite inde- 
pendent of ln(L 90 ) and O (not shown). In Fig. 4 the experimental data are 
compared with the MC- (hadron + photon) sample. The agreement for the av- 
erage log 10 (iV s /L 9 o) at fixed (—1/Rl) is not surprising because it is the result 
of the absolute calibrations, described in Section 3.2. 

MC studies have shown that the shape parameters of the lateral particle distri- 
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bution, or alternatively the age parameter in the NKG function, as determined 
in the fit to the amplitude data from the scintillation detectors, have a poten- 
tial for discriminating between photons and hadrons which may be superior to 
that of the variable 1/Rl- An approach in this direction was followed in [42]. 
In the present analysis, however, this potential could not be exploited due to 
the unsatisfactory results of the comparison between experimental data and 
MC data : no good agreement was found in the distribution of the age param- 
eter, the distribution being wider and shifted towards lower values in the MC 
data. 



4 Data Sample 



The aim of the selection criteria is to select well reconstructed showers, to 
concentrate on energies close to the threshold of the HEGRA-array and to 
retain only data taken under good observation conditions. This is achieved by 
the set of requirements (called 'standard selections' in the following): 

- ^-probability of the global AIROBICC fit : P( X 2 ) > 1% 

- zenith angle 6 < 35° 

- -0.03 in" 1 < (-1/R L ) < 0.01 in" 1 

- fitted errors of the zenith angle and of the azimuthal angle <fi : 
A0 < 0.4°, sinG • A(f) < 0.4° 

- fitted error on x and y position of the shower core : 
Ax corc < 15 m, Ay corc < 15 m 

- fitted error on ln(L 90 ) : Aln(L 90 ) < 0.3 

- ^-probability of the maximum likelihood fit to the lateral particle distri- 
bution (p s (V)) from the scintillation detectors : P(x 2 ) > 5% 

- fitted error on the total number of particles at the detector level : 
Alog 10 (i\g <0.4 

- rejection of very-high-energy showers: 

, , T s A cos0 - cos(30°) 

In L 90 < 11.0 + 0.45 ; r V-^r 

V 90; cos(10°) - cos(30°) 

This cut eliminates photon showers above ~ 320 TeV, proton showers 
above ~ 640 TeV and Fe showers above ~ 1000 TeV. The cut has to be 
applied in order to ensure consistent conditions for experimental data and 
MC data, which were generated only up to energies of 1000 TeV. 

Out of the 163 million triggered showers 19.3 million pass the standard se- 
lections. The number of triggered showers (100%) is reduced to 56%, 36% 
and 12% by requiring one after another a successful global AIROBICC fit, a 
successful fit to the lateral particle distribution p s {r) and well reconstructed 
showers. The strong reduction is explained by the large experimental errors in 
showers which are close to the energy threshold of the trigger. 
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The main characteristics of the data sample defined by the standard selections 
are (see also Figs. 1 to 4): 

- average fitted errors : 

(Ax corc ) w (Ay corc > = 4.1 m (A0) « (sin© • A0> = 0.12° 

(Aln(L 90 )> = 0.12 (Alog 10 (i\g> = 0.12 

(A(1/R L )) = 0.003 in" 1 

- average zenith angle (O) = 20° 

- average ln(L 90 ) = 9.16 (L 90 in units of no. of photons / m 2 ) 

- average log 10 (iV s ) = 4.12 

Applying the standard selections to the MC data one finds : 

- the energy region for reconstructed photon showers (10% and 90% quan- 
tiles) is 20 TeV < E 7 < 100 TeV, with an average energy of 53 TeV 

- the energy region for reconstructed hadron showers (10% and 90% quan- 
tiles) is 40 TeV < E h < 210 TeV, with an average energy of 110 TeV 

- the ratio of the number of reconstructed photon showers to the number of 
reconstructed hadron showers N 1 /N} 1 is 1/346. Note that the 7//1 flux ratio 
assumed in the MC simulation was 1/1000 at 1 TeV (see Table 1). This 
apparent suppression (by a factor of ~ 3) of hadron showers is due to the fact 
that with respect to L 90 photon showers of a given energy resemble hadron 
showers of about twice the energy, where the integrated flux is reduced by 
a factor of 2 L7 pa 3. 

The effective collection area as a function of the energy is shown in Fig. 6 for 
photon, proton and Fe showers. Defining the threshold energy as the energy 
where the effective collection area reaches 50 % of its maximum value one 
obtains for photons a threshold energy of 35 TeV, if one averages over all 
zenith angles, and 22, 28 and 45 TeV for the zenith angles 10°, 20° and 30° 
respectively. 

As can be seen from Figs. 1 to 4 the average values of the different quan- 
tities are similar for experimental data and for MC data. For MC events it 
has been verified that the average fitted error of is close to the RMS of 
(0 fitted — © generated)- Similar statements hold for other quantities. The value 
of the average fitted error (A0) ks 0.12° therefore implies that the angular 
resolution for the data sample defined by the standard selections is around 
0.12°, not including possible systematic errors. 

In the subsequent analysis the standard selections have always been applied. 
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5 Determination of Upper Limits on the Flux of Diffuse Photons 

As shown in Section 3.5 some difference between photon and hadron showers is 
seen in the plot of log 10 (iV s /L9 ) versus (—1/Rl) (Fig. 5). Under all projections 
of the 2-dimensional distribution log 10 (A r s /L 90 ) versus (—1/Rl) onto some 
axis, the one onto the v— axis, indicated by the arrow in Fig. 5, can be expected 
to yield a relative small overlap between the distributions for photons and 
hadrons. The distribution in the variable v, defined as 

t; = log 10 (AyL 90 [m- 2 ]) +23.43- (0.005 -l/R L [m]) , (1) 
will therefore be sensitive to the relative contribution of photons and hadrons. 

In the following the distributions of the variable v will be exploited for deter- 
mining upper limits of the diffuse photon background. Two approaches will 
be followed 

- Determination of a global upper limit by comparing the experimental dis- 
tribution of v with the MC expectations for primary photons and hadrons 
(Sect. 5.1). 

- Determination of an upper limit for a specific region in the galaxy by com- 
paring the experimental distributions of v in this region with that of another 
region in the galaxy (Sect. 5.2). 

5.1 Global upper limit of I^/Icr 

The distributions in v are denoted as 

Qdata = dNdata/dv for reconstructed showers of the experimental data 
g 7 = oIN^q/oIv for reconstructed MC-photon showers 
<27i = dN^c/dv for reconstructed MC-hadron showers 

Both qh and g 7 are normalized such that they correspond to the absolute 
distributions of reconstructed MC showers for an effective on-time of 1 s, 
assuming for photons and hadrons the same differential flux at 1 TeV of 0.2582 
particles / (m 2 • s • sr • TeV), which is the measured differential flux of hadronic 
cosmic rays. 

With fixed spectral indices of the photon and hadron fluxes (Table 1) the 
contribution from photons to the cosmic ray flux may be specified by the 
ratio ro = I^/Iq R of the differential photon flux and the differential hadron 
flux at some arbitrary energy E . E is chosen as 1 TeV. 
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The ratio r is then determined by fitting qdata by a superposition of q 1 and 
Qh ■ 

Qdata = T on [s] • a ■ (q h + r • q 7 ) (2) 

T on is the total on-time of the experimental data, and a and r are free pa- 
rameters to be determined in the fit. 

Using the spectral indices adopted in the MC simulation (see Table 1) r can 
be converted into a ratio r\ at any other energy E 1 (see below). If one wants to 
determine in the fit r\ instead of r , the term r -g 7 in eq. (2) has to be replaced 
by rvg 7i i with g 7i i = g 7 -r /Vi. The term qh would remain unchanged, because 
qh has been normalized to the measured differential flux of hadronic cosmic 
rays already. The choice of the reference energy Eq is therefore completely 
arbitrary. 

In case of a perfect MC simulation and an accurate measurement of the on- 
time T on the normalization factor a should come out to be 1. Leaving it free in 
the fit makes the result for tq independent of small inconsistencies in normal- 
ization between experimental data and MC data and compensates for small 
deficiencies of the MC simulation. 

The distributions g 7 , qh and qdata are displayed in Fig. 7. Due to the expected 
very small contribution from photons qdata and qh are very similar. However, 
a closer look reveals that the RMS of q da ta (0.206) is slightly lower than the 
RMS of q h (0.224). These larger fluctuations of the MC data could have their 
origin in the shower simulation or may be due to a slight overestimation of 
the measurement errors in the MC data (cf. Section 3.3.). To correct for this 
discrepancy the experimental distribution will be smeared by distributing the 
content of one bin over the same bin and neighbouring bins according to a 
Gaussian. The sigma s of the Gaussian is determined as third parameter in 
the fit. 

The experimental distribution qdata can be well fitted by a superposition of q 1 
and q h : the x 2 is 88 for 85 degrees of freedom (Fit A). In Fig. 8 the experi- 
mental distribution (points with error bars) is compared with the fitted sum 
of the contributions from hadrons and photons (dotted curve) . The fitted con- 
tribution from photons is drawn as solid curve. The distribution of the pulls 
has an average of 0.21 and an RMS of 0.99. The results for the parameters tq 
and s are 

r = 0.014 ±0.005 s = 0.080 ± 0.003 (3) 

The quoted errors correspond to 1.282 sigma, so that tq + Atq is the upper 
limit of r at the 90% confidence level, taking only the statistical errors into 
account. 
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To check the stability of the results and to assess the systematics, a and tq 
were also determined as a function of Vq, where Vq is the border line between 2 
intervals of v : v < vq and v > v . The integrated distribution of qdata, Qh an d 
q 1 in the 2 intervals were used to calculate a and r (solving two equations with 
two unknowns), keeping s fixed at 0.080. Both a and r are found to be stable 
(and consistent with the results from Fit A) in the region —0.18 < v o < 0.30. 
This is the region which contains the bulk of the showers: 74 % of qdata-, 76 
% of qh and 68 % of g 7 . For v outside of the above interval a and r are not 
stable indicating that the tails of the qdata distribution are not well represented 
by the MC simulations. 

Changing vq from —0.18 to —0.5 results in a change of r$ from 1.4% to 
— 1.5%. Setting the parameter s to zero in Fit A yields r = —3.1% with 
an unacceptable x 2 °f 465 (for 86 degrees of freedom). Using different vari- 
ables ((—l/R L ),\og w (N s /Lg )) instead of v gives either unacceptable fits (x 2 
I (no. of degrees of freedom) > 2) or r values between —3% and 0%. These 
results suggest that the tq from Fit A can be regarded as an upper limit with 
respect to systematic errors. They also imply that r$ from Fit A should not 
be interpreted as a positive photon signal, although r differs from zero by 3.6 
standard deviations. If r is fixed at zero in Fit A the x 2 value increases from 
88 to 101, corresponding to a change of the x 2 probability from 38 % to 11 %, 
which is still acceptable. The 90% confidence upper limit of r$ from Fit A is 
thus r , up i = r + Ar = 1.9 x 10~ 2 . 

The ratio r , ap i at 1 TeV is transformed to the average photon energy of the 
selected data sample of 53 TeV, using the spectral indices of the photon and 
hadron fluxes adopted in the MC simulation (see Table 1). One obtains 



The results presented so far were obtained with event samples defined by the 
standard selections. Another fit was performed using tighter cuts in and 
{-1/ Rl) : © < 15° and -0.01 m" 1 < (-1/R L ) < 0.01 m" 1 . 

The motivation for these selections is a stronger suppression of very-high- 
energy showers (0 cut) and a rejection of showers with large penetration 
depths ( (—1/Rl) cut). The latter cut preferentially rejects hadron showers 
because they exhibit larger fluctuations in the shower development. An excel- 
lent fit is obtained with x 2 — 46.8 for 56 degrees of freedom. r and r 0iUp z are 
determined as 0.012 and 0.016 respectively, yielding 



((E) = 53 TeV) < 1.4 x 10 



-2 



at 90% c.l. 



(4) 



ICR 



((E) = 31 TeV) < 1.2 x 10 



-2 



at 90% c.l. 



(5) 



ICR 
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5.2 Upper limit of I^/Icr for the inner galaxy 

The approach described in the previous section is subject to non-negligible 
systematic effects due to inconsistencies between the MC simulations and the 
experimental data. One may bypass these problems by looking for a photon 
signal using experimental data only. The MC simulation is then only needed 
to parametrize position and shape of the expected photon contribution and to 
convert the photon signal (given in terms of numbers of events) into a photon 
flux. 

Since the diffuse gamma radiation of galactic origin is expected to be con- 
centrated at the galactic disc, with a width in the galactic latitude b of a few 
degrees, it is reasonable to compare the v distribution at low \b\ (signal region) 
with that at higher |6| (background region). A larger contribution of photons 
in the signal region should then show up as an enhancement in the distribution 
for the signal region relative to the distribution for the background region. 

The acceptance of the HEGRA array depends strongly on the zenith angle © 
(due to the construction of the detectors and due to the varying overburden in 
the atmosphere) and also on the azimuthal angle (because the detector plane 
is not horizontal). In addition, because of the varying atmospheric conditions 
and because of hardware changes of the detector, the acceptance is a function 
of the time t. 

These are the reasons why a photon signal is not searched by comparing the 
v distributions in absolute scale but rather by comparing their shapes. It was 
found, however, that the shapes of the experimental v distributions depend 
on 6 and too. In principle this could be due to a non-isotropic photon 
radiation. However, the variations are very likely to a large extent caused by 
the same reasons which are also responsible for the systematic variations of 
the acceptance with 6, and t. Therefore, before comparing v distributions 
from different sky regions, which obviously are taken in different regions of 
the (Q,(j),t) space, one has to define the two v distributions such that they 
contain the same systematic effects. This is done in the following way : 

The experimental data are divided into 12 bins in cosG, 18 bins in 0, 352 
bins in t and 22 bins in v. One time bin was defined as the minimum number 
of full nights, with a total observation time of at least 5 hours, corresponding 
to a right ascension scan of > 75°. Within each bin of (cos 0,0, t) the v 
distributions of photons and hadrons are assumed to be independent of cos O, 
and t. 

With i denoting the i-th bin in the (cos 0,0, t) space and j the j-th bin of v 
the following quantities are defined : 
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no. of events from the signal region falling into the i-th bin 
of (cos Q,(j),t) and into the j-th bin of v 
total no. of signal events in the i-th bin of (cos 0, 0, t) 
normalized v distribution of events from the signal region 
in the i-th bin of (cos 6, 0, t) 
as riij but for events from the background region 
total no. of background events in the i-th bin of (cos O, 0, t) 
normalized v distribution of events from the background 
region in the i-th bin of (cos 6, 0, t) 

Because the v distribution in general depends on i, but not on cos 0, and t 
within each bin of (cos 0,0, t), information about different photon contribu- 
tions in the signal and the background region could be obtained by comparing 
the v distributions at fixed i. One may also compare v distributions summed 
over all i, provided the statistical weights of the v distributions in each bin i 
are chosen to be the same for the signal and the background region : 

The v distribution Sj and its error ASj for the signal region are defined as 

Si = E Ni ■ % = E "V (ASj) 2 = E n* (6) 

i i i 

The v distribution Bj and its error ABj for the background region are now 
constructed as 

B 3 E- v <-«-. (A^.) 2 = E^-(^) 2 (7) 

i i 

with [Awijf = rriij ■ (Mi - ra^/Mf 

By weighting the v distribution with iVj the systematic effects which are 
present in Sj are simulated in Bj. Differences between the shapes of the Sj and 
Bj distributions are now interpreted as being due to different contributions of 
photons. 

This approach is similar to the standard approach proposed by [43]. However, 
in contrast to [43], in the present analysis 

a) the background distribution is not determined from all events but only 
from those in a well defined sky region, which is different from and has no 
overlap with the signal region. Therefore signal and background distribu- 
tions (Sj and Bj respectively) are statistically independent. In addition, the 
resulting upper limit on the photon contribution is more realistic (in general 
higher) as compared to the case where signal and background regions do 



riij 

Vij Tlij j ~N i 

rriij 

Mi = Ej mj 
Wij = rriij I Mi 
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overlap. In the latter case the photon contribution would be underestimated 
in general. 

b) the background distribution Bj is not constructed by randomly generating a 
v value from the w^ distribution but rather by using for each observed event 
the full Wij distribution. In this way full use is made of the measured Wij dis- 
tribution. Both a) and b) make the error calculation simple and transparent. 

c) the sums J2i m the calculation of Sj, Bj and their errors extend only over 
those i for which M { is greater than N { . This means that in Sj and Bj only 
those (cos 0, <p, t) bins are considered for which the background distribution 
Wij is at least as well measured as the signal distribution Vij. In particular, 
if there are no background events in a certain (cos O, 0, t) bin (i.e. Mj = 0) 
then the corresponding events in the signal region (iVj) are not included in 
the calculation of Sj. The requirement M { > N { thus removes possible 
biases in the comparison of Sj and Bj. Of course, one could think of less 
restrictive conditions on the bins to be included in the sums, for example 
demanding Mj to be greater than some positive value, independent of JVj. 

The ratio r = I^/Icr a ^ 1 TeV is determined by fitting the signal distribu- 
tion Sj by a superposition of the background distribution Bj and a term Gj, 
representing a possible contribution from photons : 

S j = a-(Bj + r - C ^ (j = 1,-22) (8) 

For Gj the MC expectation g 7 (see above) for reconstructed photon show- 
ers, rebinned and normalized to J2i Ni reconstructed events, is taken, e is an 
average ratio of the number of reconstructed hadron to the number of recon- 
structed photon showers, assuming the same differential flux of photons and 
hadrons at 1 TeV. The parameters a and r are adjusted in the fit. 

The signal and background regions are defined as : 

signal region : 20° < I < 60°, and \b\ < 5° 

background region : 20° < / < 60°, and 10° < \b\ < 30° 

With this definition the average / of the events in the signal region is 48.6° and 
e is calculated as 0.348. The average I and |6| of the events in the background 
region is 45.7° and 20° respectively. The results of the fit are : 

a = 1.0007 ±0.0094 r = -(0.0003 ± 0.0030) (9) 

with a x 2 °f 26.6 for 20 degrees of freedom. The signal distribution Sj is 
compared with the fitted sum of the background distribution (a • Bj) and the 

G 

contribution from photons (a • r • — ) in Fig. 9. 
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Transforming r$ to the average energy of reconstructed photon showers from 
the signal region of 54 TeV gives the upper limit 

((E) = 54 TeV) < 2.0 x 1(T 3 at 90% c.l. (10) 

WR 

The difference (Sj — a ■ Bj) is shown in Fig. 10, together with the fitted 

G 

contribution from photons (a ■ r • — -). The distribution of the pull values 
has a mean of -0.43 and an RMS of 0.99. 



6 Discussion of the Results 



It should be noted that in the present analysis there is no way of discriminat- 
ing photon-induced showers from showers induced by electrons or positrons. 
Therefore the quoted upper limits are limits for the sum of the contributions 
from photons, electrons and positrons. Extrapolating the differential electron 
flux (using a spectral index of 3.3), as measured by [44] between 30 GeV 
and 1 TeV, and comparing with the predicted flux of diffuse photons at 50 
TeV [11] yields an e~/7 ratio in the order of 10%, at 50 TeV for the inner 
galaxy. The contribution from positrons is probably an order of magnitude 
lower than that from electrons [45]. The measured upper limits are therefore 
in good approximation upper limits of the diffuse photon flux. 

The existing measurements of Ij/Icr in the energy region from 10 GeV to 1 
PeV are displayed in Figs. 11 and 12 and compiled in Tables 2 and 3. Fig. 11 
displays the measurements which are sensitive to both the isotropic and the 
non-isotropic component of the diffuse photon flux. The measurements which 
are sensitive only to the non-isotropic component are plotted in Fig. 12. 

With one exception [29], at fixed energy all upper limits in Fig. 11 (Table 2) are 
bigger than those in Fig. 12 (Table 3). This reflects the larger systematic errors 
in measurements which are based on absolute comparisons of experimental 
data with MC predictions. In the experiment [29] the flux and the energy 
distribution of secondary photons and hadrons is measured at various depths 
in the atmosphere. It uses the different attenuation lengths for photon-induced 
showers and for the electromagnetic component of hadron-induced showers to 
determine an upper limit of I^/Icr in the energy region 5 TeV to 1 PeV. 
In [31] the experiment is critisized by claiming that the theoretical errors 
due to uncertainties in the hadronic interaction model and in the chemical 
composition may be underestimated. 

Compared to the present analysis the first measurement of a diffuse photon flux 
obtained from HEGRA-array data [31] was based on a much smaller statistics 
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of experimental data (by a factor of 250) and of MC data (by a factor of 55). In 
addition the detector consisted of 169 scintillation detectors (now 182) and 49 
Cherenkov detectors (now 97). Nevertheless the old result for r is comparable 
to the new result, which refers to slightly lower energies. In the present analysis 
more emphasis is put on a good MC tuning. Unfortunately it turns out that 
because of the increased fluctuations in the MC data (see Sect. 3. 3) and because 
of the lower energies, which imply larger errors in the measurements, the 
potential for a photon-hadron discrimination is clearly reduced. 

The upper limits in Fig. 11 (Table 2), including the result of the present 
analysis, are in general larger than the theoretical expectations by 1 to 2 
orders of magnitude. 

The measurements compiled in Fig. 12 (Table 3) are insensitive to the isotropic 
component of the diffuse photon flux because they are obtained by compar- 
ing experimental distributions in a certain sky region (usually low absolute 
galactic latitudes) with those in a background region (usually larger absolute 
galactic latitudes). 

The measurements provide a good upper limit on the non-isotropic component 
of the diffuse photon flux only if this component is negligible in the background 
region. Otherwise the upper limits will in general be too low. In the present 
analysis the background region is well separated from the signal region (5° in 
galactic latitude) in order to fullfill this condition. Note that if the IC process 
is dominating [9,13,14] a rather broad distribution of the diffuse photon flux 
in the galactic latitude is expected, in which case a good separation between 
signal and background region is even more critical. 

The upper limit from this analysis is still compatible with the theoretical pre- 
dictions shown in Fig. 12. The measurements from [46,32,55] on the other hand 
suggest that in the model by [11] the cutoff energy of the electron injection 
spectrum is well below 1000 TeV. 

The result from [32] is another HEGRA measurement in which only data from 
the scintillation detectors were used. In that analysis the lack of any gamma- 
hadron discrimination is more than compensated by the very large data sam- 
ple, collected during 2 2/3 years of data taking (day and night). In contrast 
to the data from the Cherenkov detectors, the data from the scintillation de- 
tectors depend on the atmospheric conditions (pressure and temperature) in 
a well controllable way and are practically insensitive to other conditions like 
the humidity of the air and the presence of clouds. 

It should be emphasized that the measurements compiled in Table 3 and 
plotted in Fig. 12 do not refer to the same sky regions and are therefore strictly 
speaking not directly comparable. This should also be kept in mind when 
the experimental upper limits are compared with the theoretical predictions, 
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which are for very specific sky regions. 
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Table 1 

Absolute differential fluxes = (f>o ■ (£'/TeV) _7 used for the definition of MC sam- 
ples. 



generated particle 




00 

no. of particles 
m 2 • s • sr • TeV 




spectral index 7 


sample 

is representative for 


Z 


H 


0.1091 


2.75 


H 


1 


He 


0.06808 


2.62 


He - Li 


2 - 3 





0.05182 


2.65 


Be - Si 


4 - 14 


Fe 


0.02919 


2.60 


P - Ni 


15 - 28 


H - Fe 


0.2582 




hadrons 


1 - 28 


7 


0.0002582 


2.75 


photons, e + , e~ 
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Table 2 

Measurements of I-y/IcR m the 10 GeV to 1 PeV energy range which are sensitive 
to both the isotropic and the non-isotropic component of the diffuse photon flux. 
For transforming absolute fluxes into relative fluxes a cosmic ray flux of Iqr = 
0.2582x (i?/TeV) _2 - 68 (m 2 -s-sr-TeV) _1 has been assumed. All data points are upper 
limits except those from Tien-Shan ([24]) and from EGRET [2]. The EGRET results 
are those for the isotropic component only. I and b denote the galactic longitude 
and latitude respectively. Abbreviations: ASA = air-shower array, IACT = imaging 
air Cherenkov telescope, GP1 = Galactic plane. 



experiment 


references 




Ij/ICR 


special selections 


Tien-Shan-ASA 


[24] 


> 400 TeV 


(1.0 ±0.3) x 10~ 3 


b > 50° 


emulsion 


[29] 


5 TeV - 1 PeV 


< 6 x 10~ 4 




UMC-ASA 


[46] 


> 200 TeV 


< 4.3 x 10^ 3 








> 1000 TeV 


< 4.8 x 10" 4 




Whipple-IACT 


[47] 


> 400 GeV 


< 1.1 x 10" 3 




HEGRA-ASA 


[31] 


65 - 160 TeV 


< 1.0 x 10~ 2 








80 - 200 TeV 


< 7.8 x 10~ 3 




EAS-TOP-ASA 


[48] 


> 1000 TeV 


< 7.3 x 10~ 5 








> 870 TeV 


< 1 x 10" 4 




CASA-MIA 


[49] 


> 575 TeV 


< 10" 4 




EGRET 


[2] 


10 GeV 


1.8 x 10~ 6 


out of GP1 






100 GeV 


7.0 x 10~ 6 


out of GP1 


Ooty-ASA 


[30] 


61 TeV 


< 6.7 x 10~ 2 








76 TeV 


< 3.1 x 10~ 2 








93 TeV 


< 1.9 x 10~ 2 








112 TeV 


< 7.6 x 10" 3 








175 TeV 


< 2.8 x 10- 3 








227 TeV 


< 2.0 x 10~ 3 




HEGRA-IACT 


[50] 


1 TeV 


< 8.3 x 10" 4 


GP1, 37° < I < 43°,|6| < 5° 


HEGRA-ASA 


this 


31 TeV 


< 1.2 x 10~ 2 






analysis 


53 TeV 


< 1.4 x 10~ 2 
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Table 3 

Measurements of I^/Icr m the 10 GeV to 1 PeV energy range which are sensitive to 
the non-isotropic component of the diffuse photon flux only. For transforming abso- 
lute fluxes into relative fluxes a cosmic ray flux of I C r = 0.2582 x (£'/TeV)" 2 - 68 (m 2 • 
s • sr • TeV) -1 has been assumed. All data points are upper limits except those from 
BASJE ([25]) and from EGRET [1]. The EGRET results are those for the non- 
isotropic component only. I and b denote the galactic longitude and latitude respec- 
tively. 5 and a are declination and right ascension respectively. Abbreviations: ASA 
= air-shower array, IACT = imaging air Cherenkov telescope, GPDSE = galactic 
plane drift scan experiment, GP1 = Galactic plane. 



experiment 


references 




Irs, / I ( " R 

™ 1 Oil 


special selections 


BASJE- ASA 


[251 


> 100 TeV 


8.9 x 10~ 4 


-40° < 5 < 0, 180° < a < 210° 


GPDSE 


[511 


> 900 GeV 


< 1.4 x 10~ 4 


GP1, 


|6| < 5° 






> 3 TeV 


< 1.3 x 10~ 3 


GP1, 


|6| < 5° 


BASJE- ASA 


[271 


> 180 TeV 


< 9.9 x 10~ 4 


-40° < 5 < 0, 180° < a < 210° 


UMC-ASA 


[461 


> 200 TeV 

^ AJUU J- V V 


< 8 x 10~ 5 


GP1, 


30° < I < 220°, |b| < 10° 


EAS-TOP-ASA 


[521 


> 130 TeV 


< 4 x 10~ 4 


GP1, 


\b\ < 5° 


TIBET-ASA 


[53] 


10 TeV 


< 6.1 x 10~ 4 


GP1, 


140° < I < 225°, |fr| < 5° 






10 TeV 


< 1.3 x 10~ 3 


GP1, 


20° < I < 55°, |b| < 5° 


EGRET 


[1] 


17 GeV 


2.5 x 10~ 5 


GP1, 


300° < I < 60°, |6| < 10° 






39 GeV 


3.0 x 10~ 5 


GP1, 


300° < I < 60°, |b| < 10° 


HEGRA-ASA 


[32,54] 


> 42 TeV 


< 1.6 x 10" 4 


GP1, 


0° < I < 255°, |b| < 5° 


CASA-MIA 


[55] 


140 TeV 


< 3.4 x 10~ 5 


GP1, 


50° < I < 200°, |b| < 5° 






180 TeV 


< 2.6 x 10~ 5 


GP1, 


50° < I < 200°, |b| < 5° 






310 TeV 


< 2.4 x 10~ 5 


GP1, 


50° < I < 200°, |b| < 5° 






650 TeV 


< 2.6 x 10~ 5 


GP1, 


50° < I < 200°, |b| < 5° 






1300 TeV 


< 3.5 x 10~ 5 


GP1, 


50° < I < 200°, |b| < 5° 


HEGRA-IACT 


[50] 


1 TeV 


< 2.4 x 10~ 4 


GP1, 


37° < I < 43°, |b| < 2° 


Whipple-IACT 


[56] 


> 500 GeV 


< 6.1 x 10~ 4 


GP1, 


38.5° < I < 41.5°, |b| < 2° 


HEGRA-ASA 


this 
analysis 


54 TeV 


< 2.0 x 10~ 3 


GP1, 


20° < I < 60°, |b| < 5° 
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Fig. 1. Distribution of log 10 (iV s ) for the experimental data (histogram) and for the 
MC data (open circles). 




Fig. 2. Distribution of log 10 (Lgo) for the experimental data (histogram) and for the 
MC data (open circles). 
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Fig. 3. Distribution of Ri for the experimental data (histogram) and for the MC 
data (open circles). 
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Fig. 4. Average value of log 10 (A r s /Lg ) as a function of (—1/Rl) for the experimental 
data (full circles) and for the MC data (open circles). The error bars represent the 
RMS of the \og lQ (N s / L$q) values in each bin of (—1/Rl)- 
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Fig. 5. Average value of log 10 (iV s /L9o) as a function of (—1/Rl) for photon-induced 
(full circles) and hadron-induced (open circles) showers. The error bars represent 
the RMS of the log 10 (iV s /L9o) values in each bin of (— 1/Rl). The arrow indicates 
the axis of the variable v defined in the text. 
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Fig. 6. Effective collection £1X6 £1 dbS £1 function of the energy for MC showers, after 
applying the standard selections and averaging over all zenith angles, for photon 
(solid histogram), proton (dashed histogram) and Fe showers (dotted histogram). 




Fig. 7. Distribution of the variable v for MC-photon (crosses) and MC-hadron show- 
ers (solid curve). The number on the y-axis is the number of reconstructed showers, 
after applying the standard cuts, assuming for photons and hadrons the same differ- 
ential flux at 1 TeV of 0.2582 particles / (m 2 • s • sr • TeV), the spectral indices listed 
in Table 1 and an effective on-time of 1 s. The circles represent the v distribution 
of the experimental data, multiplied by a factor 6 • 10~ 7 . 
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Fig. 8. Distribution of the variable v for the smeared experimental data (points with 
error bars). The dotted curve represents the fitted sum of the contributions from 
photons and hadrons, the solid curve the fitted contribution from photons only. 




Fig. 9. Distribution Sj of the variable v for the signal region (data points with 
horizontal bars). The polygone represents the fitted sum of the contributions from 

G 

the background region (a • Bj) and from photons (a ■ vq ■ —). 
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Fig. 10. Difference (Sj — a ■ Bj) between the v distributions of the signal and the 
background region (points with error bars). The solid curve represents the fitted 
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contribution from photons (o • ro • — -). 
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Fig. 11. Comparison of measurements of the ratio I^/Icr as a function of the energy 
in the 10 GeV to 1 PeV range. All data points are upper limits except those from 
Tien-Shan ([24]) and from EGRET ([2]). The solid line is the EGRET measurement 
of the isotropic component only. Upper limits are plotted only from those measure- 
ments which are sensitive to both the non-isotropic and the isotropic component of 
the diffuse photon background. The curves represent theoretical predictions by [11] 
with the following specifications : cutoff energy of the electron injection spectrum 
100 TeV or 1000 TeV respectively, inner galaxy (-60° < I < 60°, |6| < 10°) and 
outer galaxy (50° < I < 220°, |6| < 10°) respectively. For all predictions the spectral 
index of the electron injection spectrum was assumed to be equal to 2.0 . 
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Fig. 12. Comparison of measurements of the ratio I~/Icr as a function of the 
energy in the 10 GeV to 1 PeV range. All data points are upper limits except 
those from BASJE ([25]) and from EGRET ([1]). The solid line is the EGRET 
measurement of the non-isotropic component only. Upper limits are plotted only 
from those measurements which are sensitive to the non-isotropic component of the 
diffuse photon background only. The curves represent theoretical predictions by [11] 
with the following specifications : cutoff energy of the electron injection spectrum 
100 TeV or 1000 TeV respectively, inner galaxy (-60° < I < G0°,\b\ < 10°) and 
outer galaxy (50° < I < 220°, |fe| < 10°) respectively. For all predictions the spectral 
index of the electron injection spectrum was assumed to be equal to 2.0 . 
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